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Three independent velocities of sound can be measured along any direction of propaga- 
tion in a cubic crystal except the [100] and [111] directions. These three velocities suffice to 
determine the three elastic constants and for the [1 10] direction, the calculation of these 
constants is easy. For all other directions, the calculation is more difficult; the only existing 
method appears to be a perturbation technique developed by Neighbours. 

The present paper presents a method using exact equations and an iterative procedure 
to solve these equations and to calculate both the elastic constants and their standard 
deviations from the sound velocities and their standard deviations. The method is illus- 
trated with new data on SrTi0 3 which give Cn=3.156± 0.027, c 12 =1.027± 0.027, c 44 = 
1.215 ± 0.006X 10 12 dynes/cm 2 at 25 °C. The importance of including covariance terms in 
calculations of the standard deviations is emphasized. 

1. Introduction 

The determination of the elastic constants of single crystals from measurements of the 
velocity of sound is an extensive and active field of research and several survey papers exist 
[1, 2, 3, 4]} For any direction in a single crystal three types of sound wave may be propagated: 
one quasi-longitudinal and two quasi-transverse waves. The three corresponding velocities 
are the roots of a cubic equation, sometimes called Christoffel's equation, whose coefficients 
are complicated functions of the elastic constants and the direction cosines for the direction of 
propagation of the sound. In the case of a cubic crystal, there are only three independent 
elastic constants so that velocity measurements in a single direction suffice to completely 
determine the elastic constants provided that no two velocities are required to be equal by 
crystal symmetry. Such equality is required for the [100] and [111] directions so that measure- 
ments in one of these directions alone do not provide sufficient information to determine the 
three elastic constants. It ma} 7 happen that for some other direction two of the velocities are 
equal; the three velocities are still independent quantities, however, and measurements in such 
a direction would provide sufficient information to calculate the three elastic constants. We 
assume that the three velocities v h v 2 , v z , and their standard deviations <r h a 2 , o- 3 , have been 
measured for some direction which is specified by direction cosines /, m, n and which does not 
coincide with or closely neighbor [100] or [111]. We seek to calculate the three independent 
elastic constants c lb c 12 , c 44 and their standard deviations a n , a l2 , (r 44 . The theory leads to a 
sixth degree algebraic equation of which Cn must be a root. It may occur that more than one 
of the roots of this equation are of reasonable magnitude so that some test is needed to 
distinguish which of the roots is Cn. The direction of polarization of the quasi-transverse 
waves provides such a test and we therefore assume that the orientation of the transducer 
exciting each of the two quasi- transverse waves was also determined. 

One would usually prefer to use the [110] direction for which the calculations are easy 
and well known and for which the present method is unnecessary. However, single crystals 
of many substances are available only in very limited sizes and shapes and it may occur that 
the only available crystals do not permit measurement along [110]. Also, even if [110] is 
accessible for measurement, it may be desired to check the results by measurements in other 
directions. 



*Now at the University of Kansas, Lawrence, Kans. 

1 Italicized figures in brackets indicate the literature references at the end of this paper. 
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Neighbours and his collaborators [5, 6, 7], following a beginning by Ahrenberg [8], have 
developed an approximation method for calculating elastic constants from wave velocities 
and have applied it to the case of cubic crystals as well as several other crystal systems. In 
Neighbours' method, the equations relating the velocities and the elastic constants are expanded 
into infinite series. The first term of each equation is a simple linear combination of elastic 
constants and the first step in his self-consistent calculation is to solve for the elastic constants 
ignoring all other terms. The values so obtained are used to calculate the largest of the 
remaining terms of the infinite series and a second set of values of the elastic constants is then 
calculated considering only the elastic constants of the first term as variables. This process 
is repeated as often as necessary to obtain the desired degree of accuracy. Successive sets of 
elastic constants calculated in this manner converge to a set which satisfies the chosen finite 
portions of Neighbours' expanded form of the velocity equations. 

Presumably the next order terms could be calculated if greater accuracy were desired, 
but they apparently have not been given. The present method which uses exact equations 
thus provides a desirable alternative to Neighbours' method. The propagation-of-error theory 
could presumably be applied to Neighbours' method to calculate standard deviations, but this 
has apparently not been worked out. The present method includes such a calculation and 
this is worthy of note because the calculation involves quantities which are not statistically 
independent and serious errors can arise if the elementary propagation-of-error equations, 
which do not include co variance terms, are used. Any comparison of Neighbours' method 
with the present work should note the great utility and generality of the former which can be 
applied to crystals of any symmetry (if sufficient measurements are available) while the latter 
is strictly limited to cubic crystals. 

2. Equations for Calculating Elastic Constants 

The equations relating elastic constants to wave velocities are derived in many places; 
see, for example, Kolsky [9] for a treatment in conventional (matrix) notation and Farnell [10] 
for a brief sketch in tensor notation. The resulting secular equation can be written as an 
equation involving a 3 X 3 determinant using Farnell's notation as 

\T jk -8 Jk x\ = (1) 

where 

(2) 





X=pV z 








p = density, 








#=velocity, 
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and 

(3) 

In the last expression, the a t are the direction cosines for the direction of propagation and the 
Cmi are the elastic constants; both are referred to the crystal axes. Assuming cubic symmetry, 
writing I, m, n for the direction cosines and using the conventional matrix notation for elastic 
constants, the T Jk become 

Tn=c n l 2 +cu(rn 2 +n 2 ) (4) 

T 22 =c n m 2 +c u (l 2 +n 2 ) (5) 

Ts 3 =c n n 2 +c u (l 2 +m 2 ) (6) 

Ti2=T 21 = (c u +c 44 )lm, (7) 

ri3=r 3 i = (c 12 +c 44 )^, (8) 

and 

T 2 3 = T 32 = (c 12 +Cu)mn. (9) 
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These values can be substituted into the secular equation to give a cubic equation in x. We 
assume that for a given direction (I, m, n) the three velocities have been measured and the three 
corresponding values x u x 2 , x 3 , computed. Then for this direction the cubic equation obtained 
from the secular equation must factor into 

(x-x 1 )(x-x 2 )(x-x 3 )=0. (10) 

When this factored equation is multiplied out and the coefficients of each power of x equated 
to the coefficient of the same power in the secular equation three simultaneous equations are 
obtained. To simplify these, let 

<X = Cll, (11) 

P=cu, (12) 

y=ci2+c H , (13) 

u=x 1 +x 2 +x 3 , (14) 

v=XiX 2 +x l x 3 +x 2 x 3 , (15) 



and 

The three equations are 



and 



w=XiX 2 x z . (16) 

U = a + 2(3, (17) 

v=[d 2 +p(l 2 +m 2 )][am 2 +p(l 2 +n 2 )] 

+ [al 2 +f3(l 2 +m 2 )][an 2 +(3(l 2 +m 2 )} 

+ [am 2 +l3(l 2 +7i 2 )}[an 2 +l3(l 2 +m 2 )} 

-y 2 (l 2 m 2 +l 2 n 2 +m 2 n 2 ), (18) 

+2yH 2 m 2 n 2 -y 2 {l 2 m 2 [an 2 +f3(l 2 + 7n 2 )] 
+!V[arf+«P+n 2 )]+fflVM 2 +?(rf+n 2 )]). (19) 

The problem is to solve for a, (3, and y from a known set of values of u, v, w, I, m y and n. The 
procedure is to use eq (17) to eliminate f3 from eqs (18) and (19). Then use (18) to eliminate 
the y 2 term from (19). This leaves one equation, derived from (18), which can be solved for 
y 2 and one equation, derived from (19), which can be solved for y 3 . Cubing the first of these 
two equations and squaring the second allows the elimination of y and gives a 6th degree 
equation in a. The coefficients in this equation are very long expressions when written out in 
full and it is much more convenient to define various functions of the starting quantities (/, m, 
n, u, v, w) and so simplify the algebraic manipulations. We thus define 

k = l-l 2 , (20) 

m 1 = l-m 2 , (21) 

n! = l-n 2 , (22) 

Z 2 =3Z 2 — 1, (23) 

m 2 = 3>m 2 -l, (24) 

n 2 = 3 7l 2_ 1) (25) 

e=l 2 m 2 +l 2 n 2 +m 2 n 2 f (26) 

j = l 2 m 2 n 2 + l 2 n 2 m 2 + m 2 n% , (27) 

(j=l 2 m 2 Arl 2 n^m 2 n 2) (28) 

h=l 2 m 2 ni+l 2 n2m 1 +m 2 n 2 li, (29) 
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and 



j=l 2 m l +l 1 m2+l2ni+lin2 + m 2 ni + min2, (30) 

k=l 2 m 2 n l +l 2 n 2 m l +m 2 n 2 l h (31) 

<p=l 2 mini J rlim 2 ni J rl i min2 J (32) 

2=Zim 1 +^i + mi7ii, (33) 

r=2l 2 m 2 n 2 , (34) 

s=limin h (35) 

s 1 =l 2 m 2 n2, (36) 

9i=gJ/es u (37) 

hi = (ugk+ujj)/e—hu, (38) 

j l ^(uyk+u*qf--4vf)le--u 2 p, (39) 

ki=uk(u 2 q— 4v)/e— u s s + Sw } (40) 



p x —u 2 q— 4:V. (41) 

When (17) is used to eliminate 0, eq (18) becomes 

z?=a 2 #/4 + a^74 + ^ 2 g/4 — y 2 <?, (42) 

and eq (19) becomes 

w=A/8 + AV8 + rf#+A/8 + rT 3 -7 2 (q/+^)/2. (43) 

Substituting for 7 2 from (42) into (43) gives 

w=o? (s x e — gf) /Se + a 2 (ueh — ugk — ujf) /Se + a (u 2 ep — u 2 jk — u 2 qf+4:vf) /Se 

+ (v?es—u 3 kq-\-4:Uvk)/8e + ry d . (44) 

Computing $4r 2 y 6 /e 3 from (42) and equating to the same quantity computed from (44) gives 

a & a QJ ra 5 a 5 +a^ J ra^a 3j ra 2 a 2J raia+a =0 (45) 

where 

a =r 2 p 3 1 /e 3 —k 2 1J (46) 

a 1 =3ur 2 jp 2 Je z -2j l k 1 , (47) 

a2=Sr 2 (gpl+u 2 j 2 p 1 )/e z -2h 1 k 1 --j 2 1 , (48) 

a3=r 2 (6^> 1 +u 3 i 3 )/6 3 -2(^ 1 ^ 1 +/i 1 i 1 ), (49) 

a^&Wtf+ifaW- (2^+Aj), (50) 

a 5 =3ur 2 g 2 j/e z -2g 1 h 1 , (51) 
and 

as=rYle*-gl (52) 

The procedure for determining the elastic constants is thus as follows: Starting with the 
density, p, the velocities v i} v 2 , v 3 and the direction cosines /, ra, n first compute a?i, x 2 , # 3 from (2), 
next compute u, v, w from eqs (14) through (16), next compute the quantities defined in eqs 
(20) through (41), then compute the coefficients given by eqs (46) through (52). Using these 
coefficients plot eq (45) to determine the real, positive roots in the region of physical interest; 
if there is more than one such root choose the correct one, which is a=c n , as described below. 
Equation (45) can be plotted and the chosen root determined as accurately as desired by routine 
automatic computer techniques to save laborious hand computation. Then compute c u from 

C44=(u-eii)/2. (53) 
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Use eq (42) to compute y 2 =(c\ 2 J rCu) 2 and obtain 7. The sign of the square root is determined 
by the y 3 term in eq (43). Then compute c i2 from 

Cu=y— c u . (54) 

The only ambiguity which arises in this calculation results from (he fact that more than one 
physically reasonable root of eq (45) may occur and each such root may lead to a set of three 
elastic constants, none of which can be ruled out by the inequalities of Alers and Neighbours 
[11] or by any general physical argument. In this case, one can take one set and compute the 
polarization of each of the two quasi-transverse waves (the procedure for computing the 
polarization is given, for example, by Farnell [10]) and compare with the polarization of the 
transducer used to excite the waves. The results should agree for only one set of elastic con- 
stants so that the correct choice of root for c n can be made in this way. A second way is to 
measure velocities in a second direction in the crystal and compare the elastic constants so 
determined. Although more than one set of elastic constants may give the correct velocities 
for one direction (but not the correct polarizations) only one set should give the correct velocities 
for both directions. This second way of choosing c n avoids the need to consider the direction 
of polarization. 

3. Equations for Computing Standard Deviations of the Elastic Constants 

It is assumed that uncertainty in the density and the direction cosines may be ignored 
and that the principal uncertainty in the data is expressed by the three statistically independent 
standard deviations a h a 2) a 3 of the velocities Vi, v 2 , v s . It is easiest to divide the calculation of 
the standard deviations a n , <r 12j <iu (of c n , c i2 , Cu respectively) into two parts. First, propagation- 
of-error theory is used to calculate the variances and covariances of u, v, w. Second, these 
results are then used to calculate a n , <r 12 , and er 13 . 

The following result [12] from propagation-of-error theory is needed: Let x and y be 
statistically independent variables with known variances (variance = square of standard 
deviation). Let u and v be defined as functions of x and y and let F be defined as a function of 
u and v. Then 



where 



and 






COVI 



In some textbooks it is implicitly assumed that quantities such as u and v are statistically 
independent so that their covariance is zero and equations such as (56) and (57) with no covar- 
iance terms are given instead of the complete eq (55). The use of the complete equation is 
important in the present case; the extension of these equations from two to three variables is 
obvious. 

Application of eqs (56) through (58) gives 

v u = Vxi + 0^2 + a X3 , (59 ) 

<rl=(x 2 +X z ) 2 <r x l+ (x l + x d ) 2 a x 2 2 + (x^x^a^, (60) 

oi=(x 2 X z ) 2 cr x l+ (x 1 x 3 ) 2 a x 2 2 + (xiX 2 ) 2 a x l, (61) 

Cov (u, v) = (x 2 +x 3 )a x l+(x 1 +x s )a x l+(x 1 + x 2 )a x 2 d , (62) 

Cov (u, w)=x 2 x Z (j x \-\-XiX 3 (T X2 J rXiX 2 (T xz , (63) 
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and 

Gov (v, w) = (x 2 +x 3 )x 2 x z <r x 2 1 + (x 1 +x d )x 1 x 3 a x 2 2 + (Xi+X 2 )X 1 X 2 (T Z 1 J (64) 

where 

(J x i = 2 pVidiy <T x2 = 2pV 2 (T 2 , <JxZ= z 2pV z d Z . 

To use eq (55) for the standard deviations of the elastic constants one must compute the 
partial derivatives of these constants with respect to u, v, w. These partials of c n can be 
obtained by taking derivatives of eq (45). The resulting expressions involve partials of the 
coefficients a t and these in turn involve partials of the quantities defined in eqs (38) through (41). 
Let subscripts u, v, w denote partial differentiation with respect to u, v, w respectively. Then 
from eqs (38) through (41) 

h lu =(gk+jf)/e-h, (65) 

jiu=2u[(jk + qf)/e-p), (66) 

k lu =k(3u 2 q-4:v)/e-3u 2 s, (67) 

Piu = 2uq, (68) 

j lv =-4f/e, (69) 

k lv =—4uk/e, (70) 

Pi.= -4 (71) 

k lw =S (72) 

and the following are all zero: h lv , h lw , j lw , and p lw . We shall next require the partial derivatives 
of the a t which are given by differentiating eqs (46) through (52) with* the result for the u 
derivatives 

a 0u =3r 2 p 2 1 p lu /e*—2k 1 k lu , 

aiu = 3r 2 (jp 2 1 + 2ujp 1 p lu )/e d —2j 1 k lu —2k 1 j lu , 

a 2u = 3r 2 (2gp 1 p lu +2uj 2 p l +u 2 j 2 p lu )/e 3 —2h 1 k lu —2k 1 h 1 u-2j 1 j lu , 

a2u=r 2 (6gjp 1 + 6ugjp lu +3u 2 f)/e*—2(g 1 k lu +h^ lu ^^ 

a4 U =3r 2 (2ugj 2 +g 2 p lu )/e 3 —2g l j lu —2h 1 h lu , 

a 5u = 3r 2 g 2 j/e 3 -2g 1 h lu , 

and with a du equal to zero. For the v derivatives the result is 

a 0v = — 12r 2 pl/e 3 — 2k x k lV) 
a lv =—24r 2 ujp l le*—2j l Jc lv —2k 1 j lv , 
a 2v ^-12r 2 (2gp 1 +u 2 f)/e"-2h 1 k lv -2j 1 j lv , 
a 3v =—24r 2 ugj/e z —2(g l k lv +h 1 j lv ) } 
a Av = — 12r 2 g 2 /e s —2g 1 j lv , 

and with a 5v and a &v equal to zero. For the w derivatives the result is 

a 2w = — 16h 1 , 
a Zw = — l6gi } 

and with a 4tn a bw , and a 6w equal to zero. Define 

D=QaQCi 1 + 5a 5 c 1 \+4ta 4 c^ J r3a d c 2 l + 2a 2 c n + a l (88) 
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(76) 
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(79) 
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(80) 
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(84) 




(85) 




(86) 




(87) 


V 



then differentiation of eq (45) gives 

Cnu= — (a 5u c; ) l +a iu c^+as u c 1 \+a 2 uCn + a lu Cn+ao u )ID } (89) 

Cnv= — ((i4vC l \ + as v Ci l + a 2v c l 2 1 +a lv c n + aov) /D, (90) 

and 

Cnw= — (a3wCn+a 2w Ci\-\-ai w c n +a 0w )/D, (91) 

where we have assumed ZMO. For the [110] direction D=0 and a different treatment, described 
below, is required. The derivatives of c 4A are obtained from (53) and are 

Cuu=(l-c nu )l2, (92) 

c U v=— CuJ2, (93) 

and 

Cuw=— Cuw/2. (94) 

The derivatives of c u are obtained from eq (54) and so involve derivatives of y. From (42) we 
have 

y 2 = ((jc 2 l +ujc n + p 1 )4e, (95) 

so that 

(y 2 )u=(2c Al ()c nu +jCn+ujCnu+Piu)/4:e, (96) 

(y 2 )v=(2c n gc llv +ujc lu -4:)/4:e, (97) 

and 

(y 2 )w=(2c n g+uj)c Uw /4e. (98) 

Now use y u =(y 2 )u/2y and eq (54) to obtain 

Ci2u=(y 2 ) u /2y-c iiu , (99) 

cv2*=(y 2 )v/2y-c i4v , (100) 

and 

c^w=(y 2 )w/2y-c iiw . (101) 

The expression for o- n is then 

o'i 2 i=Cii^or^+Cii^+Ci 1 ^o-^ + 2ci lM c 1 i, cov (w, v)+2c nu c Uw cov (u, w)+2Cn D Ci ltt , cov (p, w). (102) 

The equations for a M and <r 12 are the same with the subscripts on c n changed to 44 and 12 
respectively. 

The procedure for obtaining the desired standard deviations is thus straightforward 
although tedious. One begins by calculating the variances and covariances of u, v, w from eqs 
(59) through (64). Then compute in succession the quantities given by eqs (65) through (102). 

Following eq (91) we noted that D=0 for the [110] direction. This can be seen as follows: 
For the [110] direction r is zero by eq (34) and the t 3 term drops out of eq (43). Then eq (45) 
simply consists of the square of all the terms in (43) except ry*. Let 

B=a\/8 + a 2 uh/8 + au 2 p/8 + uh/8-y\af+uk)/2-w. (103) 

Then eq (45) for the [110] direction is B 2 =0, and D is then D=2B -^, and therefore D = 0. 

The foregoing statistical treatment thus fails for any direction for which c n is a double 
root of eq (45). This appears to be true only for the [110] direction (we have already noted 
that the [100] and [111] directions are not suitable for the method of this paper) but the writers 
have not been able to construct a proof. 
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For the [110] direction the following results are well known and easily obtained from the 
treatment of Kolsky [9], for example. If x x corresponds to the longitudinal wave, x 2 to the 
transverse wave with displacement parallel to [001], and x 3 to the transverse wave with dis- 
placement parallel to [110], then 

Cu=Xi — x 2 +x z , (104) 

Ci2=x 1 — x 2 —x z , (105) 

and 

c u =x 2 . (106) 

These give 

(7 ll ==o 'l2 = o 'zi + 0':r2+ o 'z3> (107) 

and 

a u = (Tx2- (108) 

We have assumed throughout this paper that errors in thickness and density can be 
neglected in comparison with errors in transit time t. If standard deviations were assigned to 
the thickness and density as well, the quantities x t would not be statistically independent and 
two modifications of the foregoing treatment would be required. First, eqs (56) through (58) 
would have to be used with the thickness, density, and transit times as independent variables 
to give the variances and covariances of u, v, w. The calculations using eqs (65) through (102) 
would then go as before. Second, eqs (107) and (108) would have to be replaced by appropriate 
expressions in terms of the standard deviations of the thickness, density, and transit times 
derived from (56) through (58) and (104) through (106). No attempt has been made to allow 
for uncertainties in the orientation; such allowance should probably be made but appears to 
be an exceedingly difficult problem. 

4. Procedure for Computing the Elastic Compliances and Their Standard 

Deviations 

The foregoing results form a complete unit, giving the elastic constants, c ih and their 
standard deviations. The present section may be skipped unless it is desired to put the results 
in terms of the elastic compliances, s^. The calculation of the 8 ti from the c tJ is trivial but the 
calculation of the standard deviations of the s fj is more difficult and serious errors can result 
if the covariance terms are not taken into account. There appears to be no discussion of this 
problem in the literature on elastic constants, so we present the method for the cubic case. 

The well-known equations for the elastic compliances of a cubic crystal in terms of the 
elastic constants are 

*n=(cii+c u )/c, (109) 

Si 2 =— cu/c, (110) 
and 

s Ut = 1/cu, (111) 
where 

c=(c 11 -c 12 )(cn+2c 12 ). (112) 

To compute the standard deviations of the s tj one can apply eq (55) which requires the 
covariances of the c ijt To compute these covariances one might try to use eq (58) which would 
be wrong because x and y were assumed to be statistically independent. We require the more 
general formula 

fJPn . dFdG 2 ,dFdG 2 ,/dFdG,dFdG\ , , , 11Q . 

cov (F, G)=^ 55 «l+ Tv ^ d+(^ ^+^ ^) cov (u, v) (113) 

for the covariance [12] of F and G which are defined in terms of quantities u and v which are not 
statistically independent. Writing cov (11, 12) for the covariance of c n and c 12 we have 
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COV (11, 12)=C nu C 12u ^l + CnvC 1 2v(Tl + Cn w Ci2w^l+(CnuCi2v + CnvCi2u) COV (u, V) 

+ (c nu C 12w + C U wCi2u) COV (U, W) + (c U vCi2w + Cu w Ci2v) cov (v, w) (114) 

The expressions for cov (1 1, 44) and cov (12, 44) are identical except for the appropriate changes 
of subscripts, but we shall not need to calculate these latter two covariances. Writing Sn,u for 

^ — and similarly for other partials we have 

a, u = g - fel+ ^ 1+ H (115) 

g „, 12 == c + (Cll+Cl2) ^ 2+Cll) . (116) 

Ci 2 (2cii+Ci 2 ) / 117 x 

£12,11=-- - ~2 — — > U 1 '; 

(Cn + Ci2)(4c 12 — Cu) niQ N 

§12,12 = 2 ' (Ho) 

8« M = -\ldU, (119) 

and with §11,44, §12,44, §44,11, §44.12 all zero. Letting o-§u represent the standard deviation of §n 
and similarly for the other s ih 

we have 

^ 11 =§ii,no'ii+§ii,i2^i2+2§ii,ii§ii,i2 cov (11, 12) (120) 

0-? ]2 =§12, 11^11 + ^12, 12^12 + 2§12,11§12. 12 cov (11, 12), (121) 

and 

<r su =cjcl. (122) 

For the [110] direction the covariances of c u and c 12 is best calculated directly from (104) and 
(105) using (58). The result is 

cov (11, l2) = a 2 xl +<T 2 z2 +al 3 . (123) 

Thus, the procedure for calculating the standard deviation of the s t j for any direction, 
including the [100] and [1 1 1], is to use eqs (120) through (122) evaluating the partial derivations 
from eqs (115) through (119). The situation considered in the present paper, using only 
information obtainable from measurements in a single direction, compels us to exclude [100] 
and [111] from the present considerations. For [110] cov (11, 12) is given by (123); for any 
other direction, it is given by (114). 



5. Measurements on SrTi0 3 

The writers carried out a series of measurements on a boule of strontium titanate, kindly 
supplied by the National Lead Company, to check the method. The density value [18] of 5.116 
g/cm 3 was used. Flats were first ground to give the maximum thickness between parallel 
faces permitted by the shape of the crystal. A series of measurements was taken and analyzed. 
The crystal was then recut normal to the [110] direction and a second series of measurements 
was then made. The measurements were all made with 10 Mc/s X-cut and AC-cut quartz 
crystals 0.25 in. in diameter. A commercial pulsed oscillator was used to drive these transducers. 
The echoes were observed on a dual trace oscilloscope simultaneously with a 1 Mc/s standard 
frequency signal. The results are sumarized in table 1. 
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Table 1. Data on SrTi0 3 at 25° C 



Parameter 


Symbol 


First 
direction 


Second direc- 
tion, [110] 


Direction 


I 

m 
n 


0. 67905 

. 27160 

. 68200 

5. 3045 

13. 15=b0. 10 

22. 08=k. 05 

23. 15±. 07 

3. 330±. 051 

1. 181±. 005 

1. 074=b. 006 


0. 70711 


Cosines 


. 70711 


Length... ... _ cm.. 


. ooooo 

4. 5670 


Transit times, 10 -6 sec 


h 

X\ 
X2 
X3 


11. 36±0. 05 


Xi, 10 12 dyne/cm 2 


18. 71±. 07 
20. 02±. 09 
3. 308=b. 029 




1. 219±. 009 
1.064±.01O 



6, Results 

The data are summarized in table 1. The calculation of c n by the method of eqs (20) 
through (52) was programmed for an automatic computer using an iterative procedure of solving 
eq (45) which gives the real roots to four places in the interval 0.0 to 5.0 X10 12 dynes/cm 2 . 
For the first direction, there is a single root at a=3.162 and a single root at a = 3.497X10 12 
dynes/cm 2 . For the second direction there is a double root at a=3.153 and a double root at 
a = 3.462 X10 12 dynes/cm 2 . A complete set of elastic constants was computed for each of these 
possible choices of c n and the results are compared in table 2. The choice of 3.497 and 3.462 
leads to a disagreement in c l2 which is outside the experimental error. The choice of 3.162 
and 3.153 gives consistent sets of constants. The latter choice is also known to be correct 
because 3.153 is obtained from eq (104) when x 2 and x 3 are properly distinguished by the polari- 
zation of the corresponding sound waves. 

For the first direction, the standard deviations of the c tj were determined by the method 
of eqs (59) through (102) using an automatic computer. Equations (107) and (108) were used 
for the second direction. The final values for the c tj were computed by averaging the results 
for the two directions weighted by the reciprocals of the squares of the standard deviations. 
The Sij values were then computed from eqs (109) through (112) and their standard deviations 
from eqs (113) through (123). 

Table 2. Comparison of possible sets of elastic constants 



Constants 


First direction 


Second direction 


Computed from correct choice of root for cn 


Cn 

Ci2-- -- - -- - _- - 


3. 162±0. 052 
1.035=1= .052 
1.212=1= .007 


3. 153±0. 032 
1.024± .032 




1. 219=1= .009 






Computed from incorrect choice of root for cn 


Cn 

Cl2 

C44 


3.497 
1.200 
1.044 


3. 462=1=0. 032 
1.024=1= .032 
1. 064=1= .010 



The c n and c 12 values of Bell and Rupprecht [14] agree within experimental error with 
the results of the present work as shown in table 3; the c x2 value of Poindexter and Giardini 
[15] also agrees within experimental error, but their c n value differs from the present result by 
much more than twice the standard deviation and so is significantly different from our result. 
The other workers' c 44 values lie on either side by slightly more than twice the standard deviation 
of our value, but are probably within the combined experimental error of their determination 
and ours. The writers feel that the c n values of Poindexter and Giardini should be rejected 
and that the remaining data show reasonable agreement. 

The anisotropy of a cubic crystal depends on the quantity 8 defined by 



5=2s n —2s 12 —s u . 
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(124) 



The Young's modulus, Y fj and shear modulus, 6 fy are given as a function of the usual spherical 
polar angles by 

_ = s n __5 sin 2 cos 2 6- (6/4) sin 4 6 sin 2 20, (125) 

1 r 



and 



_= 544 -l.25 sin 2 69 cos 2 (9+ (5/2) sin 4 6 sin 2 20. 



<*/ 



(126) 



These quantities are plotted in figure 1 which shows that SrTi0 3 comes close to being isotropic; 
Young's modulus varies by only 1 percent and the shear modulus by 5 percent. 

Table 3. Comparison with elastic constants of SrTiOz at 25° C determined by other workers 



Constants 


Poindexter 

and 
Giardini a 


Bell and 
Rupprecht *> 


Present 
work b 


Percent 

difference, 

last two 

columns 


Cii 

Ci2— - -_._____ 


3.48 
1.01 
1. 19 

3.3 

-0. 74 
8.4 


3. 181 
1.025 
1.236 

3.729 

-0. 909 

8.091 


3.156±0.027 
1.027± .027 
1.215± .006 

3. 772± . 023 

-0.926± .010 

8.233± .040 


0.79 
0.19 


Cu -_ ._ 


1.73 


8n~ --- 


1.14 


Sl2 . 


1.84 


Su 


1.73 







All dj in units of 10 12 dyne/en^; sa in units of 10— 13 cm2/dyne. 

a Sij determined by resonance met hod and cy by matrix inversion. 

b en determined by pulse velocity method and 8%j by matrix inversion. 



Figure 1. Young's modulus, 

Y f , and the shear modulus, 

Gf, as a function of orien- 
tation. 

The colatitude, 0, is the angle between 
the [001] direction and the direction of 
measurement. The azimuth, <f>, is the 
angle between the [100] direction and 
the projection of the direction of 
measurement on the (001) plane. The 
subscript / indicates that the elastic 
moduli are for a free specimen which is 
under no constraint. 
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7. Summary 

(1) Velocity measurements in a single direction in a cubic crystal provide enough informa- 
tion to determine the three elastic constants, c ih except for the [100] and [111] directions which 
are therefore excluded from consideration in this paper. 

(2) For the [110] direction the computation of the elastic constants and their standard 
deviations is simple and is given in eqs (104) through (108). 

(3) For all other directions the calculations are much more complex. The general pro- 
cedure for the elastic constants (applicable also to the [110] direction) is given in eqs (2) through 
(54). The general procedure for the standard deviations (not applicable to the [110] direction) 
is given in eqs (59) through (102). 

(4) The procedure for calculating the elastic compliances, s ih and their standard devia- 
tions is given in eqs (109) through (123). Throughout the statistical treatment the covariance 
terms are included and their importance is emphasized. 

(5) The method is applied to SrTi0 3 and results in good agreement with previous workers 
are obtained. 



The writers thank W. E. Tefft of the National Bureau of Standards for bringing eqs (55) 
and (58) to their attention, and M. D. Beals of the National Lead Company for supplying the 
SrTi0 3 specimen. 
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Selected Abstracts 

A new type of computable inductor, C. H. Page, ./. Research 

NBS MB (Math, and Math. Phys.), No. 1 (Jan.-Mar. 1968), 
75 cents. 

The mutual inductance analog of the generalized Thompson- 
Lamparo theorem (for cross-capacitances) is developed. 
An infinitely long cage of five parallel wires can yield an 
absolute inductance of 



henries per meter. End-effects of order l/l 2 occur in a 
finite cage, but can be reduced to order l// 4 by using eight 
wires. 

The eight wire cage has the advantage of overdetermined 
relations among the inductances to be measured, allowing 
an estimate of experimental error in the calibration of a 
standard. Errors due to faulty cage geometry are shown 
to be of the order of 1 in 10 7 . 

Numerical computation of the temporal development of 
currents in a gas discharge tube, \V. Borsch-Supan and IT. 
Oser, J. Research NBS 67B (Math, and Math. Phys.), No. 1 
(Jan.-Mar. 1968), 75 cents 

The behavior of electrical currents in a gas discharge tube 
including space charge effects is investigated by numerical 
integration of the governing nonlinear partial differential 
equations. Both stationary solutions and the temporal 
development, under the influence of space charge effects, 
are considered. The numerical integration is done by a 
difference method. Il is found that the truncation error can 
be greatly reduced by comparison with formal solutions for 
constant fields. The discussion is essentially restricted to 
the more mathematical questions. 

Equipment for single crystal growth from aqueous solution, 
J. L. Torgesen, A. T. Horton, and C. P. Saylor, ./. Research 
NBS QIC (Eng. and Instr.), No. 1 (Jan.-Mar. 1968), 75 cents. 
The growing of large single crystals of high quality from 
solution requires the precise control of supersaturation and 
the avoidance of thermal amd mechanical shock. Uniform 
growth conditions and cleanliness need to be maintained. 
Good seed crystals are necessary and the accidental introduc- 
tion or generation of new nuclei should be prevented insofar 
as possible. In the apparatus and equipment assemblies 
here described, the crystal-growth bath is designed for 
uniform growth conditions and the exclusion of contamina- 
tion. The support for the crystals provides for convenient 
mounting of the seeds, holds the crystals firmly, and allows 
their easy removal with minimum damage. A new reversing 
rotation mechanism promotes equal washing of the crystal 
surfaces and achieves uniform temperature and composition 
of the solution without, at the same time, producing significant 
mechanical stresses in the crystals. The temperature 
controller gives regulation of the temperature an order of 
magnitude more sensitive than those hitherto used and pro- 
vides for stepless change of temperature. The crystals are 
thus free from liquid inclusions found to result from sudden 
acceleration of growth. The electrolytic conductance of 
ionic solutions may be used as a precise and sensitive measure 
of solution concentration and supersaturation. The crystal- 
growth procedures which are reported have resulted in the 
production of very good single crystals of a considerable 
variety of chemical phases. 

An automatic multichannel correlator, R. F. Brown, Jr.# 
/. Research NBS 67C (Eng. and Instr.), No. 1 (Jan.-Mar 
1963), 75 cents. 

Very early in our infrasonic Research Program it become 
apparent that the only way to differentitate between local 



meteorological pressure variations and the infrasonic signals 
from distant sources was through the relative time (^ arrival 
of the pressure disturbance at pickups spaced several miles 
apart. The "Multichannel Correlator" is an analog computer 
that was designed to detect these infrasonic signals. 
At each pickup the pressure variations are converted to FM 
audio signals which are transmitted over leased telephone 
lines to the correlator. At the correlator these FM signals 
are demodulated and recorded on a }{ inch per minute multi- 
channel magnetic tape. The recorded data frequency 
range is 0.02 to 1.0 c/s and the pressure range is 0.1 to 50 dynes 
per square centimeter. 10 minute blocks of data are read 
from the tape and amplitude equalized to a constant rms 
value. A programmed time delay device produces the 
equivalent delays appropriate to a continuous horizontal 
azimuth search for velocities in the range of 280 to 400 meters 
per second. The correlator output is a continuous record of 
the average of the rectified sum of the delayed channel 
signals. 

Elastic constant porosity relations for polycrystalline thoria, 

S. Spinner, V. P. Knudsen, and L. Stone, ./. Research NBS 
67C (Eng. and Instr.), No. 1 (Jan.-Mar. 1963), 75 cents. 
The relations between the elastic constants and porosity for 
about 300 thoria specimens have been determined. Both 
Young's and the shear modulus for each specimen were ob- 
tained by a dynamic resonance method. From these two 
moduli, Poisson's ratio was computed. The tendency for 
ail I hese elastic constants to decrease wit h increasing porosity 
is in qualitative agreement with theory. However, the 
decrease in the experimental values was greater than would 
be expected from the theory. This greater decrease for the 
experimental values is attributed to the fact that the speci- 
mens do not conform to the idealized assumptions of the 
theory. 

New fast-opening, large-aperture shutter for high-speed 
photography, E. C. Cassidy and I). II. Tsai, -/. Research 
NBSWC (Eng. and Instr.),' No. 1 (Jan.-Mar. 1968), 75 cents. 
This shutter consists essentially of a metallic foil in a capacitor 
discharge circuit. The opening action is obtained when the 
foil is compressed by the electromagnetic forces which ac- 
company the heavy surge current through the circuit, during 
a transient discharge. A shutter made up of two foils in a 
loop arrangement may be opened to an area 1 in. x 3 in. in 
less than 45 jusee. The factors affecting the design and 
operation of this shutter are analyzed. These factors include 
the circuit parameters, the energy input to the foil, and the 
materials and the size of the foil. Some experimental results 
are also given. 

Perturbation method in a problem of waveguide theory, 

D. Fox and W. Magnus, J. Research NBSQTD (Radio Prop.), 
No. 2 (Mar. -Apr. 1963), 70 cents. 

The reflection coefficient for the basic mode in a widening, 
straight, two-dimensional waveguide is computed for small 
wave numbers by using the perturbation method with the 
electrostatic case as the unperturbed case. The problem is 
treated as a perturbed infinite system of inhomogeneous 
linear equations, and it is shown that the matrix of the un- 
perturbed system (which corresponds to the electrostatic 
case) can be inverted explicitly by using conformal mappings 
and physically unrealistic modes. Questions of convergence 
are discussed, and other examples for application of the 
method are indicated. 

Some wave functions and potential functions pertaining to 
spherically stratified media, C. T. Tai, /. Research NBS 

67D (Radio Prop.), No. 2 (Mar.-Apr. 1963), 70 cents. 
The wave functions pertaining to a bilinearly stratified 
medium are presented in this paper. Solutions to the equa- 
tion V • [«(r)vV]=0 have been investigated for several 
profiles of *(/). An analysis is given to the equation 
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V • [" • v^]=— p/e which arises from the formulation of 
the quasi-static electric field in a homogeneous anisotropic 
medium. 

Radiation from a plasma-clad axially-slotted cylinder, 

W. V. T. Rusch, J. Research NBS 67D (Radio Prop.), No. 2 
(Mar -Apr. 1963), 70 cents. 

Expressions have been developed for the complete radiation 
pattern from a finite axial slot on an infinite circular cylinder 
covered with a homogeneous, uniform plasma sheath. The 
resulting field expressions are functions of both the longitude 
angle, <£, and the polar angle, 0. Numerical results have 
been obtained for a thin lossless plasma to determine the 
variation of the field pattern with plasma frequency, plasma 
thickness, slot width, and cylinder radius. These results are 
presented graphically to illustrate various properties of the 
radiating system. 

Report of the 47th National Conference on Weights and 
Measures, 1962, NBS Misc. Publ. 244 (November 23, 1962), 
75 cents. 

A report of the proceedings of the forty-seventh National 
Conference on Weights and Measures, held in Washington, 
B.C., June 4, 5, 6, 7, and 8, 1962, and attended by State, 
county, and city weights and measures officials. 167 p. 

Radiation quantities and units, International Commission on 
Radiological Units and Measurements (ICRU) Report 10a, 
1962, NBS Handb. 84 (Nov. 14, 1962), 20 cents. 
This Handbook presents definitions of twenty-three funda- 
mental radiation quantities and units. It resulted from a 
three-year study by the Ad Hoc Committee on Quantities 
and Units of the ICRU. It includes new names for certain 
quantities and clarified definitions for others. It presents a 
system of concepts and a set of definitions which is internally 
consistent and yet of sufficient generality to cover present 
requirements and such future requirements as can be foreseen. 

Controlled temperature oil baths for saturated standard cells, 

P. H. Lowrie, Jr., NBS Tech. Note 141 (August 1962), 25 cents. 
Two oil baths for the temperature control of saturated stand- 
ard cells have been designed and fabricated at the Boulder 
Laboratories of the National Bureau of Standards for opera- 
tion at 28° C and 35° C respectively. Short term control to 
better than ±0.001° C with day-to-day variations no greater 
than 0.002° C has been achieved with the use of a mercury- 
toluene thermoregulator incorporating a temperature antici- 
pating device. The circulating system limits temperature 
gradients in the oil to less than 0.001° C across any 10 inch 
section. The baths incorporate pre-heat and drain tanks as 
well as the main temperature regulated tank to facilitate the 
insertion and removal of cells and to minimize oil spillage. 

Applications of resistance thermometers to calorimetry, 

G. T. Furukawa, Book, Temperature, Its Measurement and 
Control in Science and Industry III, pt. 2, 317-328 (1962). 
The importance of the resistance thermometer in the accurate 
measurement of both temperature and the heat leak of the 
calorimeter is discussed. The final accuracy of the determi- 
nation of heat capacity is shown to be dependent upon the 
accurate and consistent measurement of heat input to the 
sample and the corresponding rise in temperature. The 
various heat-capacity calorimeters used in the range from 10 
to 400 °K are briefly described with emphasis upon the 
applications of resistance thermometers, the methods for 
calibrating them, and the problems associated with the 
design of calorimeter vessels. Comparison is made of the 
thermometric properties of platinum, copper, indium, lead 
and gold-silver alloy. The need for high relative accuracy 
in the measurement of AR is emphasized. The various 
temperature scales used in calorimetry are compared and 
their applications are described. 

Absorption spectrum of carbon vapor in solid argon at 4° 

and 20 °K, R. L. Barger and H. P. Broida, J. Chem. Phys. 
37, No. 5, 1152-1153 (Sept. 1, 1962). 

Absorption spectra have been observed photoelectrically for 
C 3 and C 2 in solid argon at temperatures of 4° K and 20° K. 
A molecular beam of carbon from a Knudsen cell at 2550° K 



was deposited simultaneously with argon onto a quartz 
surface, cooled either with liquid helium or with liquid hydro- 
gen. The spectra consist of 21 lines and a continuum be- 
tween 3797A and 4221A attributed to C 3 , and two lines at 
5208 A and 5269 A attributed to C 2 . Measured wavelengths, 
half-widths, peak absorption intensities and /-values are 
given. An /-value of 2X 10 -4 was found for the 5269 A line 
of the C 2 and the 4102A line of C 3 . 

Vibrational-rotational spectroscopy, H. C. Allen, Jr. and 
W. B. Olson, Ann. Rev. Phys. Chem. 13, 221-240 (1962). 
The literature on molecular spectroscopy has increased 
significantly each year with the increased availability of 
high-quality commercial infrared spectrometers. This year 
is no exception. It is quite impossible to give a complete 
review of all the work in the space available, hence the authors 
will take advantage of their prerogative of discussing a 
limited number of topics. We shall, thus, stress work of a 
general nature on theoretical advances, fine structure of 
simple molecules, vibrational spectra, intensities of gases, 
and computer applications. For the spectra of specific 
molecules the reader is referred to Chemical Abstracts. 

Realistic diatomic potential function, H. W. Woolley, /. 
Chem. Phys. 37, No. 6, 1307-1316 (Sept. 1962). 
An examination of recently published data for diatomic 
molecules as obtained by the Rydberg-Klein-Rees method 
has suggested a useful form for the potential energy function. 
This term provides an inverse power dependence at large 
distance and an extrapolation to infinity at the extreme of 
close approach. A three parameter potential obtained by 
fitting in lowest approximation is compared with empirical 
data on co e X e and a e . Relations are also given to permit an 
arbitrarily close fit of observed spectroscopic constants. 
Theoretical formulas covering seven new Yu's of the Dunham- 
Sandeman series for the rotating-vibrator are included. 
Potential curves for H 2 and HF are shown as examples, 
including effects of previously determined higher order 
spectroscopic constants. Other applications of the new 
potential are also discussed. 

Optical constants of aluminum in vacuum ultraviolet, R. 

La Villa and H. Mendlowitz, Phys. Rev. Letters 9, No. 4, 
149-150 (Aug. 15, 1962). 

Values of the frequency dependent complex dielectric constant 
e(co) of aluminum in the photon energy region of 12-16 ev 
were derived from the characteristic electron energy loss 
spectra. A two-parameter Drude model with N, the number 
of "free" electrons per atom, as 2.6 and r, the relaxation time 
as 1.1X 10~ 15 sec. was found to give a good approximation 
to the optical properties of aluminum. 

Microwave spectrum of methyldifluoroarsine, L. J. Nugent 

and C. D. Cornwell, /. Chem. Phys. 37, No. 3, 523-534 
(Aug. 1962). 

The spectrum of CH 3 AsF 2 was investigated over the range 
10-33 KMc/sec. Three Q-branches (K=5-4, 6-5, 7-6) 
and two # transitions (J=2- 3, K = 2-3; /=3-4, if = 3-4) 
were fitted by a Hamiltonian which included terms for over- 
all rotation, internal rotation, nuclear quadrupole coupling, 
and centrifugal distortion. The reciprocal moments of 
inertia are A = 5414.49± 0.04 Mc/sec, £=5381. 28± 0.03 
Mc/sec, and C=3871.88±0.03 Mc/sec. The 6-axis is 
perpendicular to the plane of symmetry. The height of the 
internal rotation barrier (assumed to be sinusoidal) is 466 ± 10 
cm -1 . With the bond length rcH= 1.095 A assumed, the 
spectral data yield r AsC = 1.930± 0.010 A, r A «,p= 1.734T 0.005 
A ZHCH=109 o 9'±l°, ZCAsF=95°40'±14', ZFAsF=96°18' 
±16'. The quadrupole coupling constants for As 75 are 
eQq C c=- 220± 4 Mc/sec and Vo = (q bb - q aa )lq C c = - 1.20± 0.05. 
The indicated uncertainties of all parameters are limits of 
error obtained directly from the estimated maximum im- 
precision in the measured frequencies. 

New wavemeter for millimeter wavelengths, R. W. Zimmerer, 
Rev. Sci. Instr. 33, No. 8, 858-859 (Aug. 1962). 
A new wavemeter of simple design is described. The principle 
of operation makes use of a new development in physical 
optics. The actual performance of the devise was measured 
and compared with the theory. 
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Electron energy losses in solids and thei 
electron diffraction diagram, L. Marlon, ./ 
17, Suppl. B-II, 68-73 (1961). 
The electron diffraction diagram is conside 
of view of the influence of characteristic e 
electrons in the solid. Some of the presei 
losses are reviewed, together with the expt 
leading to a new viewpoint. While the 
characteristic losses on the average electr 
gram is negligible, for precise measuremei 
the losses and a corresponding correction 
or experimental) is needed. 
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Optical detection of microwave transitions in electronically 
excited CN produced by a chemical reaction, R. L. Barger, 
H. P. Broida, A. J. Estin, and H. E. Radford, Phys. Rev. 
Letters 9, No. 8, 345-346 (Oct. 15, 1962). 

Population inversions and non-equilibrium populations pro- 
duced by chemical reactions should be useful for molecular 
structure studies by double resonance and possibly for maser 
or laser action. Atomic flames of nitrogen atoms mixed 
with hydrocarbons provide a particularly good example of a 
chemical reaction which produces population inversion. A 
TEon microwave cavity of volume 40 cm 3 , from 8750 to 
10,200 Mc and having a Q of approximately 2X 10 4 , was used 
with a power input of \{ to 1 watt. At each of the frequencies 
9855, 9740, and 8810 Mc, a resonant increase in the optical 
emission was observed; the increases had ralative magnitudes 
of 5, 4, and 1, respectively. 

Sum rules for vibrational-rotational energy levels including 
centrifugal distortion, H. C. Allen, Jr., and W. B. Olson 
J. ehem. Phys. 37, No. 2, 212-214 (July 15, 1962). 
Sum rules are developed for the vibrational-rotational energy 
levels of an asymmetric rotor. These sum rules contain the 
corrections for centrifugal distortion correct for terms of the 
order P 4 in the Hamiltonian. If the 1 M terms of the Hamil- 
ton ian are neglected these sum rules reduce to those derived 
for a rigid asymmetric rotor by Mecke. 

Rotational perturbations in CN I. Zero-field theory, optical 
Zeeman effect, and microwave transition probabilities, ft. E. 

Radford and II. P. Broida, Phys. Rev. 128, No. 1, 231-242 
(Oct. 1963). 

The theory of rotational perturbations in doublet states of 
diatomic molecules is developed and is applied to the long- 
standing problem of perturbations in the violet bands of the 
CN radical. Analysis of the available optical data yields 
experimental values of two reduced matrix elements: 



(n|AL y |2)= 
(n|BL y |2) = 



=-0.39±0.04 cm- 
=0.011 ±0.001 cm- 



where the uncertainties shown are estimated upper limits 
These account satisfactorily for all the observed line shifts in 
the violet spectrum. These two parameters are then used to 
predict the Zeeman effect of the rotational perturbations, 
and these predictions are checked against new observations 
of the CN violet specturm at magnetic field strengths up to 
28 kgauss. Finally, the zero-field analysis is used to examine 
the feasibility of a proposed microwave resonance investiga- 
tion of the perturbed excited states of CN. 

The system of B12O3B2O3, E. M. Levin and C. L. McDanieL 
,/. Am. Ceram. Soc. 45, No. 8, 355-360 (Aug. 1962). 
The phase diagram for the system Bi 2 03-B 2 0.3 has been 
determined experimentally. The melting point of Bi 2 3 has 
been redetermined as 825° C with an estimated overall 
uncertainty of about ±3° C, and the molal heat of fusion of 
Bi 2 3 , calculated from the slope of the liquidus curve, is 2050 
cal per mole. The system contains a body-centered cubic 
phase of approximate composition 12Bi 2 3 -B 2 03, which melts 
incongruently at 632° C. Four congurently melting com- 
pounds exist iu the system: 2Bi 2 3 B 2 03, 3Bi 2 3 -5B 2 3 , 
Bi 2 03-3B)O;j, and Hi >() :r I B2O3, with melting points, respec- 
tively, of 675°, 722°, 708°, and 715° C. The Bi 2 3 -4B 2 3 
compound exhibits a sluggish transformation at 696° C. 
Compositions containing up to 97.5 wt% (85 mole %) 
Bi 2 3 can be partly or totally quenched to glass. Indices of 



the quenched glasses are greater than 1.74. A region of 
liquid immiscibility extends at 709° C from almost pure 
B2O3 to 19.0 mole % Bi 2 3 . The extent of immiscibility 
theoretically calculated agrees with the experimentally cle- 
termined value when 1.20 A is used for the ionic radius of 
Bi 3 +. 

Wavelengths, energy levels, and pressure shifts in mercury 
198, Y. Kaufman, /. Opt. Soc. Am. 52, No. 8, 866-870 (Aug. 
1962). 

The vacuum wavelengths of 27 lines in the spectrum of Hg 198 
have been measured relative to the international standard of 
length, the Kr s6 line at 6057.8021 1 A by photographic Fabry- 
Perot interferometry. These measurements were made using 
Hg 198 electrodeless lamps containing argon gas at pressures of 
}i, 3 and 10 mm Hg. Energy level values have been derived 
for each of the Hg 198 sources. On the basis of these values, 
the energy level shifts per unit pressure of argon have been 
calculated. The suitability of the Hg 198 electrodeless lamp 
as a source of wavelength standards for interferometric 
measurement of length and wavelength is discussed. 
The standard source used was a Kr 8fl hot cathode lamp 
operated according to the conditions proscribed by the 
International Conference on Weights and Measure's in 1960. 

A network transfer theorem, G. F. Montgomery, IRE Trans. 
Audio. AU-10, No. 3, 88 (May-June 1962). 
For a linear, passive, reciprocal two-port network, the 
forward, open-circuit voltage transfer ratio is equal to the 
reverse, short-circuit current transfer ratio. 

Acoustical interferometer employed as an instrument for 
measuring low absolute temperatures, ./. Acoust, Soc. Am. 
34, No. 8, 1145-1146 (Aug. 1962). 

Values of absolute temperatures at 2 and 20° K have been 
determined from experimental measurements of the speed of 
sound as a function of pressure in helium gas. The acoustical 
interferometer was the instrument employed in the measure- 
ments, and the accuracy achieved in the experiment indicates 
that sonic thermometry at low tempertures may be competi- 
tive with other conventional thermometry techniques. 

Kinetics of Cs 1 ' desorption from tungsten, M. D. Scheer and 
J. Fine, J. Chem. Phys. 37, No. 1, 107-113 (July 1962). 
The mean absorption lifetime (r 3 ) of Cs + on tungsten has 
been measured in the 1 000-1 200 °K temperature range, using 
a pulsed beam technique. Under conditions of low surface 
coverage and with either Cs or Csl as beam materials, 
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was obtained. The heat of desorption can be calculated as 
the energy required to remove an isolated Cs + ion from the 
surface of an electrical conductor. The presence of an 
absorbed contaminating layer, arising from residual vacuum 
gases, decreased the Cs + — W binding energy by 0.5 ev and 
increased the pre-exponential factor by about two orders of 
magnitude. Anomolous results were obtained when CsCl 
was used for a beam material, suggesting a reaction between 
the surface tungsten atoms and atomic chlorine. 

Structure and structure imperfections of solid (3 -oxygen, 

E. M. Horl, Acta, Cryst. 15, No. 9, 845-850 (Sept. 1962). 
Electron diffraction studies were carried out on thin films of 
solid oxygen. The structure of 0-oxygen was investigated 
and a rhombohedral structure of the space group P3i21 
found. Its corresponding hexagonal cell has the dimensions 
a= 3.307 ±.008A, c=11.256± .015A and contains three 
molecules parallel to the hexagonal axis at the positions: 
(i, i, 0), (-h 0, |) and (0, -J, f). Faults in the stacking 
sequence of the (OO.l)-layers were observed and the experi- 
mental results compared with the Paterson treatment of 
growth faults in fee crystals. 

Neutral meson photoproduction from complex nuclei, II. A. 

Schrack, J. E. Leiss, and S. Penner, Phys. Rev. 127, No. 5, 
1772-1783 (Sept. 1, 1962). 



207 



The angular distributions of neutral mesons produced by 
170 Mev bremsstrahlung on carbon, aluminum, copper, 
cadmium, and lead have been obtained through the coincident 
detection of the decay photons by a scintillation counter 
system. The experimental data have been analyzed by 
means of a Monte Carlo prediction based on an impulse- 
approximation elastic coherent production model. The 
nuclear density distributions which were used in the synthesis 
of those predictions which were in best agreement with the 
experimental data from this experiment have been compared 
to the charge density distributions inferred from electron 
scattering experiments. The values of the rms radii of 
nuclear matter obtained in this experiment are, within the 
limits of error of the experiment (about 2.X 10~ 14 cm), con- 
sistent with the values of charge distribution radii obtained 
in electron scattering. 

The spin-independent part of the neutral meson photopro- 
duction cross section used in the synthesized predictions is 
obtained using the dispersion theory of Chew, Low, Gold- 
berger, and Nambu. The dependence of the cross section on 
the meson-nucleon phase shifts used is investigated and a set 
of phase shifts is presented that results in a correct prediction 
of the absolute cross sections observed in this experiment; 
however, this particular set of phase shifts is not unique. 

Rotation-vibration spectrum of matrix-isolated hydrogen 
chloride, L. J. Schoen and D. E. Mann, /. Chem. Phys. 37, 
No. 5, 1146-1147 {Sept. 1, 1962). 

An infrared vibration-rotation spectrum of hydrogen chloride 
molecules isolated in solid argon matrices at 1:500 dilution 
has been observed at 4 and 20 °K in the region of the funda- 
mental. From the contraction in spacing of the P(l) and 
R(0) it has been deduced that the HC1 molecules undergo 
quantized but hindered rotation in cavities having cylindrical 
symmetry. 

Other NBS Publications 

Journal of Research 67B (Math, and Math. Phys.), No. 1 
(Jan-Mar. 1963) 75 cents. 

Evaluation of a generalized elliptic-type integral. L. F. 

Epstein and J. H. Hubbell. 
An algorithm for obtaining an orthogonal set of individual 

degrees of freedom for error. J. M. Cameron. 
Recognition of completely mixed games. A. J. Goldman. 
A new type of computable inductor. C. H. Page. (See 

above abstract.; 
Numerical computation of the temporal development of 

currents in a gas discharge tube. W. Borsch-Supan and 

H. Oser. (See above abstract.) 
Tables of genera of groups of linear fractional transformations. 

H. Fell, M. Newman, and E. Ordman. 

Journal of Research 67C (Eng. and Instr.), No. 1 (Jan.- 
Mar. 1963), 75 cents. 

Method for calibrating a standard volt box. B. L. Dunfee. 
Stability of residual thiosulfate in processed microfilm. C.I. 

Pope. 
Equipment for single crystal growth from aqueous solution. 

J. L. Torgesen, A. T. Horton, and C. P. Saylor. (See 

above abstract.; 
An automatic multichannel correlator. R. F. Brown, Jr. 

(See above abstract.) 
Elastic constant — porosity relations for polycrystalline thoria. 

S. Spinner, F. P. Knudsen, and L. Stone. (See above 

abstract.) 
An oxygen partial pressure warning instrument. L. 

Greenspan. 
New fast-opening, large-aperture shutter for high-speed pho- 
tography. E. C. Cassidy and D. H. Tsai. (See above 

abstract.) 
Equations for the radiofrequency magnetic permeameter. 

C. A. Hoer and A. L. Rasmussen. 

Journal of Research 67D (Radio Prop.), No. 2 (Mar.-Apr. 
1963), 70 cents. 

The protection of frequencies for radio astronomy. R. L. 
Smith- Rose. 



Radar reflections from the moon at 425 Mc/s. G. H. Millman 
and F. L. Rose. 

Sunset and sunrise in the ionosphere: effects on the propaga- 
tion of longwaves. J. Rieker. 

Correction of atmospheric refraction errors in radio height 
finding. W. B. Sweezy and B. R. Bean. 

Empirical determination of total atmospheric refraction at 
centimeter wavelengths by radiometric means. A. C. 
Anway. 

Propagation of radiofrequency electromagnetic fields in geo- 
logical conductors. V. Fritsch, translated from German 
by A. P. Barsis. 

WWV reception in the arctic during ionospheric disturbances. 
G. E. Hill and J. R. Herman. 

Height-gain for VLF radio waves. J. R. Wait and K. P. 
Spies. 

Perturbation method in a problem of waveguide theory. 
I). Fox and W. Magnus. (See above abstract.) 

Some wave functions and potential functions pertaining to 
spherically stratified media. C. T. Tai. (See above ab- 
stract.) 

Radiation from a plasma-clad axially-slotted cylinder. 
W. V. T. Rusch. (See above abstract.) 

Two- and three-loop superdirective receiving antennas. 
B. W. Seeley. 

Hallen's method in the problem of a cavity-backed rectangu- 
lar slot antenna. J. Galejs. 

Relative convergence of the solution of a doubly infinite set 
of equations. R. Mittra. 



Periodicals received in the Library of the National Bureau of 

Standards, July 1962, N. J. Hopper, NBS Mono. 57 (Nov. 

23, 1962), 25 cents (Supersedes NBS Circ. 563 and the 1st 

supplement to NBS Circ. 563). 
Corrosion of steel pilings in soils, M. Romanoff, NBS Mono. 

58 (Oct. 24, 1962), 20 cents. 
Hvdraulic research in the United States 1962, H. K. Mid- 

dleton, NBS Misc. Publ. 245 (Oct, 26, 1962), $1.00. 
Quarterly radio noise data, March, April, May 1962 and cor- 
rigendum for Technical Notes 18-1 through 18-11, W. Q. 

Crichlow, R. T. Disney, and M. A. Jenkins, NBS Tech. 

Note 18-14 (Aug. 9, 1962), 50 cents. 
Mean electron density variations of the quiet ionosphere, 

No. 8— October 1959, J. W. Wright, L. R. Wescott, and 

D. J. Brown, NBS Tech. Note 40-8, (Sept. 1962) 35 cents. 
Synoptic radio metrology, B. R. Bean, J. D. Horn, and L. P. 

Riggs, NBS Tech. Note 98, (Oct. 1962) 50 cents. 
Bibliography on direction finding and related ionospheric 

propagation topics 1955-1961, O. D. Remmler, NBS Tech. 

Note 127, (Oct. 1962) 60 cents. 
Equatorial spread F, W. Calvert, NBS Tech. Note 145 (Aug. 

1, 1962) 60 cents. 
The thermodynamic properties of helium from 6 to 540° R 

between 10" and 1500 psia, D. B. Mann, NBS Tech. Note 

154A (Jan. 1962), 50 cents. 
The energy parameter B for strong blast waves, I). L. Jones, 

NBS Tech. Note 155, (July 1962) 25 cents. 
Emission stabilization of thermionic diode noise sources, 

M. W. Randall and M. G. Arthur, NBS Tech. Note 160 

(Sept. 1962), 15 cents. 
Evaluation of unexpectedly large radiation exposures by 

means of photographic film, W. L. McLaughlin, NBS 

Tech. Note 161, (Aug. 1962) 15 cents. 
Spectrophotometric determination of hydroperoxide in 

diethyle ether, W. C. Wolfe, Anal. Chem. 34, No. 10, 

1328-1330 (Sept. 1962). 
Exchange behavior of kaolins of varying degrees of crystal- 

linity, W. C. Ormsby, J. M. Shartsis, and K. H. W 7 oodside, 

J. Am. Ceram. Soc. 45, No. 8, 361-366 (Aug. 1962). 
Analysis of the hvdroxyl radical vibration rotation spectrum 

between 3900 A and 15000 A, A. M. Bass and D. Garvin, 

J. Mol. Spectry. 9, No. 2, 114-123 (Aug. 1962). 
Microwave measurements in the NBS electronic calibration 

center, R. E. Larson, Inst. Elec. Engrs. 109, pt. B, suppl. 

No. 23, 644-650 (1962). 
Effect of additives on silver iodide particles exposed to light, 

G. Burley and D. W. Herrin, J. Appl. Meteorol. 1, No. 3, 

355-356 (Oct. 1962). 
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A modulated sub-carrier technique of measuring microwave 

attenuation, G. E. Schafer and R. R. Bowman, Inst. Elec. 

Engrs. 109, pt. B, suppl. 23, 783-786 (1962). 
A survey of microwave power-measurement techniques 

employed at the National Bureau of Standards, G. F. 

Engen, Inst. Elec. Engrs. 109, pt, B, suppl. No. 23, 734- 

739 (1962). 
Symposium on teaching of materials: Aspects of material 

behavior significant to engineers, J. M. Frank land, J. Eng, 

Mech. Div.-Proc. Am. Soc. Civil Engr. Paper No. 3246. 

75-81 (Aug. 1962). 
Excitation of modes at very low frequency in the earth- 
ionosphere wave guide, J. R. Wait, J. Geophvs. Res. 67, 

No. 10, 3823-3828 (Sept. 1962). 
The growth and roots of electronic computers, J. L. Little, 

Student Quart, and EE Digest, Joint Publ. of Am. Inst. 

Elec. Engr. and IRE 1, No. 1, 64-70 (Sept. 1962). 
Acoustical interferometer employed as an instrument for 

measuring low absolute temperatures, G. Cataland and 

H. H. Plumb, J. Acoust., Soc. Am. 34, No. 8, 1145-1146 

(Aug. 1962). 
The Hall effect, D. O. Webster, Student Quart. & EE Dig. 

p. 82-86 (Sept. 1962). 
Dimensional changes in complete dentures on drying, wetting, 

and heating in water, G. C. Paffenbarger and W. T. 

Sweeney, J. Am. Dental Assoc. 65, No. 4, 495-505 (Oct. 

1962). 
A variable-parameter direct-current switching filter, G. F. 

Montgomery, Proc. IRE 50, No. 9, 1986 (Sept. 1962). 
The vapor phase radiolysis of propane-^/ 8 in the presence of 

other hvdrocarbons, S. G. Lias and P. Ausloos, J. Chem. 

Phys. 37, No. 4, 877-883 (Aug. 15, 1962). 
Synthesis of an immittance function with two negative im- 
pedance converters, S. B. Geller, IRE Trans. Circuit 

Theory CT-9, No. 3, 291 (Sept. 1962). 
A note on the propagation of electromagnetic pulses over the 

earth's surface, J. R. Wait, Can. J. Phys. 40, 1264-1269 

(1962). 
Surface behavior of silver single crystals in fused sodium 

chloride, J. Kruger and K. H. Stern, J. Elect rochem. Soc. 

109, No. 10, 889-894 (Oct. 1962). 
The spontaneous martensitic transformations in 18% Cr, 

8% Ni steels, R. P. Reed, Acta Met. 10, 865-877 (Sept. 

1962). 
A high speed pyrometer, G. A. Hornbeck, Book, Tempera- 
ture, Its Measurement and Control in Science and Industry, 

III, pt. 2, 425-428 (Reinhold Publ. Corp., New York, 

N.Y., 1962). 
The viscous heating correction for fiscometer flows, E. A. 

Kearsley, Trans. Soc. Rheology VI, 253-261 (1962). 
Strengthening of hot work die steels, C. R. Irish and S. J. 

Rosenberg, Trans. Quart. 55, No. 3, 613-623 (Sept. 1962). 
Correlation of factors influencing the pressures generated in 

multi-anvil devices, J. C. Houck and U. O. Hutton, Am. 

Soc. Mech. Engrs. Paper 62-WA-254 (1962). 
The temperature dependence of flow and fracture charac- 
teristics of an age-hardenable alloy, W. D. Jenkins and 

W. A. Willard, Trans. ASM 55, No. 3, 580-598 (Sept. 

1962). 
The effect of experimental variable including the martensitic 

transformation on the low-temperature mechanical prop- 



erties of austenitic stainless steels, C. J. Guntner and R. P. 

Reed, Trans. ASM 55, No. 3, 399-419 (Sept. 1962). 
A new ionospheric multipath reduction factor (MRF), R. K. 

Salaman, IRE Trans. Commun. Systems CS 10, 220-222 

(June 1962). 
The effect of electrolytes on I lie sodium chromotropism of 

hi\-(meso-2,3-diaminobutane)-nickel (II) ions, I). L. 

Leussing, J. Harris, and P. Wood, J. Phys. Chem. 6<>, 

1544-1546 (Aug. 1962). 
Characteristics of resistance strain gages, R. L. Bloss, Book' 
Semiconductor and Conventional Strain Gages, ed. 

Mills Dean, III and R. D. Douglas, ch. VII, 123-142 

(Academic Press, Inc., New York, N.Y., Oct. 1962). 
The photodisintegration of bismuth and the lend isotopes, 

E. G. Fuller and E. Hayward, Nuclear Phys. 33, 431-448 

(1962). 
Comparison of the relative acidic strengths of the isomeric 

dinitrophenols in benzene and water, M. M. Davis, J. 

Am. Ceram. Soc. 84, 3623-3627 (1962). 
The basis of the measurement system, A. V. Astin, Proc. 

IRE 50„ No. 5, 614-616 (May 1962). 
Predicting the performance of long distance tropospheric com- 
munication circuits, A. P. Barsis, K. A. Norton, and P. L. 

Rice, IRE Trans. Commun. Systems CS-10, 2-22 (Mar. 

1962). 
A neuron model which performs analog functions, D. R. 

Boyle, (Proc. San Diego Symp. Biomedical Engr., San 

Diego, Calif., R. 256-263 .lime 19-21, 1962). 
Scattered radiation from large Cs n7 sources, L. Costrell, 

Health Phys. 8, No. 5, 491-498 (Oct. 1962). 
Standard tests for electrical properties, A. H. Scott, SPE J. 

18, 1375-1378 (Nov. 1962). 
The industrial significance of National Bureau of Standards 

research, G. M. Kline, Mod. Plastics 40, No. 2, 149 (Oct. 

1962). 
Comments on "A new precision low level bolometer bridge," 

G. F. Engen, Proc. IRE 50, No. 9, 1997 (Sept. 1902). ' 
Solar particles in interplanetary space, C. S. Warwick, Skv 

and Telescope XXIX, No. 3, 133-136 (Sept. 1962). 
The kinetics of the heat precipitation of collagen, J. M. 

Cassel, L. Mandelkern, and D. E. Roberts, J. Am. Leather 

Chem. Assoc, LVII, No. 11, 556-575 (Nov. 1962). 
Chromium plating by thermal decomposition of dicumene 

chromium, W. H. Metzger, Jr., Plating 49, No. 11, 1176 

(Nov. 1962). 
Oxidation of high polymers, M. Tryon and L. A. Wall, 

Book, Autoxication and Antioxidants II, ch. 19, p. 919 

(Interscience Publ. Inc., New York, N.Y., 1962). 
The use of a Venturi tube as a quality meter, R. V. Smith, 

P. C. Wergin, J. F. Ferguson, and R. B. Jacobs, J. Basic 

Eng. 84, 411-412 (Sept, 1962). 
Dislocation loops in deformed copper, A. W. Ruff, Jr., 

Fifth Intern. Congress for Electron Microscopy, p. J-10 

(Academic Press, Inc., New York, N.Y., 1962). 

^Publications for which a price is indicated are available 
fro?n the Superintendent of Documents, U.S. Government 
Printing Office, Washington 25, D.C. (foreign postage one-fourth 
additional). Reprints from outside journals and the NBS 
Journal of Research may often be obtained directly from the 
authors. 
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